An Index Theorem for the Majorana Zero Modes in Chiral P-Wave Superconductors 
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We show that the Majorana fermion zero modes in the cores of odd winding number vortices of a 
2D px + ipy-paired superconductor is due to an index theorem. This theorem is analogous to that 
proven by Jackiw and Rebbi for the existence of locahzed Dirac fermion zero modes on the mass 
domain walls of a ID Dirac theory. The important difference is that, in our case, the theorem is 
proven for a two component fermion theory where the first and second components are related by 
parity reversal and hermitian conjugation. 
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The odd winding number vortices in a spinless (spin- 
polarized) px + ipj,-paired superconductor or superfluid 
(for brevity we shall use "superconductor" to denote both 
in the rest of the paper) trap zero-energy bound states 
in the cores. The second-quantized operators creating 
these zero energy excitations are the self-hermitian, Ma- 
jorana fermion operators. This property endows the vor- 
tices with non-Abelian statistics [1-3], which means that 
under braiding of any two vortices the total wavefunc- 
tion transforms as a vector in a finite dimensional Hilbert 
space. Because of the non-Abelian statistics, this type of 
superconductor has been proposed to support topological 
quantum computation [2, 4, 5]. For example it has been 
shown that the PfafFian = f fractional quantum Hall 
state [1, 6] is a spinless + ipy superconductor of the 
"composite fermions" [1] . The braiding of non-Abelian 
quasiparticles in such a state[7] can, in principle, be ex- 
ploited to build a quantum computer that is immune to 
environmental errors [8] . 

The root of the non-Abelian statistics is the presence of 
Majorana fermion zero modes in the vortex cores. Such 
zero modes were first proposed to exist in the vortex 
cores of the px + ipy-paired superfluid A-phase of He- 
lium 3 in Ref. 9. Traditionally this peculiar kind of vor- 
tex core state was demonstrated by explicitly solving the 
Bogoliubov-de-Gennes (BdG) equations[l, 5, 10, 11]. In- 
terestingly, the zero-energy solution is found only in the 
case when the winding number of the vortex is an odd in- 
teger (for even winding number vortices, there is no zero 
energy mode). 

The zero mode in the vortex core discussed above is 
very analogous to the zero energy domain wall states 
of polyacetylene[14]. In that case, by writing down an 
ansatz for the dimerization-order -parameter profile, one 
can also demonstrate the existence of the zero energy 
domain wall solution by explicitly solving the mean-field 
equations[15]. Quite satisfactorily, this domain wall state 
was shown to be a condensed matter realization of the 
zero mode associated with the mass solitons of a ID Dirac 
theory investigated earlier by Jackiw and Rebbi[12, 13]. 



The Jackiw and Rebbi soliton solution is a simple exam- 
ple of an index theorem where fermionic zero modes can 
be used to count the topological defects (or magnetic flux 
quanta) of a background order parameter (or magnetic 
field). 

In this paper, we ask what the analogous index the- 
orem is for the vortex Majorana fermion zero modes of 
a Px + ipy superconductor. We proceed by mapping the 
2D vortex problem on an effective ID problem by per- 
forming angular momentum decomposition with respect 
to the center of the vortex. In this way, we can show that 
for odd winding number vortices there exists a unique an- 
gular momentum channel in which the following Hamil- 
tonian describes the low energy quasiparticle excitations: 



Hm = I dx 



(1) 



In Eq. (1), m{~x) — ~m{x) is a spatially varying mass 
term that changes sign at x = (the location of the 
domain wall), and x\x) is a two component field given 
by x^(a;) — {c\x), c(— a;)) , with c{x) being a spinless 
fermion field. Note the important difference of Eq. (1) 
with the Dirac theory. 



Hd = I dx 



ivpip^ Uzdxtp + m{x)il'^ Ox^p , (2) 



where %Ij^{x) = [f\(x), f\{x)) with fi^2{x) being two in- 
dependent fermion fields. It is because of such difference, 
the zero energy quasiparticles localized on the mass do- 
main walls of Eq. (1) are Majorana fermions, while those 
localized on the domain walls of Eq. (2) are ordinary 
fermions. As we will show later, for even winding number 
vortices, such angular momentum channel does not exist. 

We first briefiy review the derivation of the Jackiw- 
Rebbi zero mode for Eq. (2). The quasiparticle operator 

= [ dx [<ly^{x)fl{x)+cj,2{x)fl{x)] 



2 



satisfies [H^q^ = eq^. That implies the foUowing Dirac 
equation for the two component wavefunction (j)^{x) = 
{4>i{x),(j)2ix)): 



—ivpazdx4>{x) + Uxm{x)4>{x) = e(j){x) 



(3) 



First we note that because CTy anticommutes with ax and 
CTz, if (f){x) is an eigenfunction with eigenvalue e, C7y(j){x) 
is also an eigenfunction with eigenvalue — e. As a result, 
the e = solutions of Eq. (3) can be made a simultaneous 
eigenstate of ay. Let (j)o{x) denote such a solution and 
ay4>o{x) = X(j)o{x). Set e = and left-multiplying Eq. (3) 
by iaz we obtain 



dx(t>o{x) = —m{x)(l)o{x), 

Vp 



which implies 



(4) 



For m{x) = ±sign(a;)|m(x)|, Eq. (4) is normalizable for 
A = ^1. In this way we have proven that for each sign 
change of m{x) there is a single zero energy mode. 

Now let us consider a uniform (i.e. with no vortices) 
2D px + ipy superconductor. The fermionic mean-field 
Hamiltonian is given hy K + Hqp , where 

k 

Hop = -Ao Y^{kx + iky)clc\ + h.c. (5) 

k 

Here K is the kinetic energy term, Hqp is the pairing 
term and = fc^/2m — ep with ep the Fermi energy. 

With the purpose of treating the electronic state of a 
single vortex in mind, we change to a new representation 
where the fermion operators are expanded in angular mo- 
mentum channels, 



Ck 



1 



E 



with m an integer. The commutation relation {ck, c|,} = 
5^{k — p) implies 

{Cmk, 4p} = 6m,nS{k - p). 

Inserting them in Eq. (5), doing the 6^ integral, and lin- 
earizing around ^f, we get. 



K = 



1 



(27r) 



.A 

/ dq (vpq) c\^„Cm,q, 
J-A 



(6) 



where q = k — kp, vp is the Fermi velocity, and A is a 

momentum cut-off. So the kinetic energy term separates 
into uncoupled angular momentum channels indexed by 



the integer angular momentum m. Using similar manip- 
ulations to decompose the pairing term, we find, using 
^_k = TT + 6u and kx + iky = |fc|e'^''. 



OP 



27r2 



..A 

dqcos{miT)cl^gc\_^g+h.c. (7) 



Putting Eq. (6) and Eq. (7) together we have, for each 
pair of m and 1 — m, the following massive Dirac theory 



Hri 



-I 



dx 



- iv'ptjjl^azdxipm + moipl^axtlJr, 



In the above, ijj}n{x) 

y^rt 



IS 



the Fourier transform of 

^m,qy^i-m,,q), Vp — anO TOq — gi^ . 

Consequently, all fermionic quasiparticle excitations are 
gapped. 

Next we consider the fermionic Hamiltonian for a single 
winding number vortex located at tlie origin. The kinetic 
energy part of the Hamiltonian remains the same as in 
Eq. (6). Let us now consider the pairing term. In order to 
describe the spatial dependence of the superconducting 
order parameter, we start with the real space description. 



H 



IP 



..Jd^Rj 



d're^''^h{R)g{r)cl^^cl_, + h.c.. 

(8) 

Here, R and r are the center-of-mass and the rela- 
tive coordinates of the Cooper pair, respectively. h{R) 
and are the amplitude and the phase of the su- 
perconducting order parameter, and ^(r) is the Fourier 
transform of (px + iPy)- In the vortex core, h{R) 
(1 — e~~) with ^ the coherence length. Substituting 
cjj^ir = 27r5^k4e'''-*^'^'"^ in Eq. (8), we end up with 
two spatial integrals, ^(k — p) = J d'^rg{Y)e'^^^~^^''' = 
{kx - Px) + i{ky - Py) and 



/(k + P) = J (i'i?e*^^/i(i?)e'(^+P)-^. 



(9) 



In order to evaluate /(k-|-p), we first note that, I{Re{k+ 

p)) = e*^/(k + p, where Rg is the operator that rotates 
(k -|- p) by an angle 6 in the momentum space. It follows 
that, /(k-l-p) = e'^''+p7(|k-hp|). To evaluate /(|k-|-p|) 
we choose (k + p) along the y-axis. Performing the 0r 
integral which produces — 27ri J_i(|k-|-p|i?), where J_i is 
the Bessel function of the first kind of order —1 [16], and 
then performing the R integral which produces |^]f^p|2 x 
we find. 



i^lP = -(27r)3^AoX^^^±i^ 



{Px +iPyf 



k,p 



k + p|' 



h.c. 



Finally, using angular momentum expansion of the 

fermion operators and noting that a function of |k + p| is 
periodic in (0k — ^'p) and hence can be Fourier expanded 



3 



as 



|k + p|3 

we rewrite Hip as, 

ifiP = -(27r)2zAoEk,pE„.™„ 



■P J 



(10) 

Performing the 6'k and 9p integrals in Eq. (10) for the 
two different terms, we get tlie conditions 

mi = 2 — m, TO2 = m, and mi ~ — m, m2 ~ 2 + m 

respectively. Since |^,^p|a depends on cos(6'k; — 9p), and 
hence is even in ((?k — ^^p), u„i{k,p) = U-m{k,p). In 
addition, since |^,_^^p|3 is symmetric in the k p ex- 
change, Um{k,p) = Um{p,k). Using these, we perform 
the transformation m — > — m and interchange k and p in 
the second term of Eq. (10) to arrive at the simple form 
for the pairing term, 

iIiP = -2iAo^ / dkdpk'^^/kpum{k,p)cl_^ ,^4^ p+h.c. 

(11) 



In the above pairing Hamiltonian, each angular mo- 
mentum channel m is coupled to channel (2 — m). The 
only exception is the channel to = 1 which is decoupled 
from the rest. In this channel, writing k = k^ + q and 
p = kp + q' for small q and q' , and noting that the coeffi- 
cient of c| ^, must be odd in (g — q') from the fermion 
anticommutation relation (otherwise, the integrals over q 
and q' in the pairing term below will give zero since the 
fermion bilinear changes sign under q q'), the Hamil- 
tonian for the TO = 1 channel takes the form, 

^1 = T^~yl / '^^ + iAo / dqdq' A{q - q') 

(12) 

where A(q — q') is an odd function of {q — q'), which, in 
leading order, is simply proportional to {q — q')- Defining 
a two component fermion operator xj = (ci^jCi,,), 
Hi becomes Eq. (1) where xH^) ~ (c^(^)i c(— a;)) is 
the Fourier transform of Xg ^^'^ m{x) is the Fourier 
transform of iAoA{q). Due to the odd nature of A{q), 
the mass term satisfies m{—x) = —m{x). 

Now we prove the index theorem for Eq. (1). By writ- 
ing the quasiparticle operator as 



7 



j dx[rii{x)x\{x) + ri2{x)xi^2{x)] 



where x}i{x) = c^(x) and xJC^) = '^(^■'^) demanding 
that [-ff, 7^] = €7^, we obtain the following equation for 
r{^{x) = {r]i{x),r]2{x)), 



-ivF<Jzdxri{x) + m{x)axr]{x) = -?7(a;) 



(13) 



The rest of the proof for the zero mode is identical to that 
in the Dirac fermion case. To see that the zero energy 
quasiparticle is a Majorana fermion we recall that the 
zero energy solution of Eq. (13) is an eigenstate of ay. 
Explicitly it takes the form 



r]{x) = e^p 



t/4 



\/2 



Here we have chosen a particular global phase for the 
(7y eigenvector. Like the Dirac fermion zero modes, the 
value of A here is ±1 when the mass profile satisfies 
m{x) = =Fsign(x)|TO(a;)|. The corresponding quasipar- 
ticle operator takes the form. 



pj7r/4 



r 



1 J dxe^^Srniv)<iy ^_lZ(ci{x)-ic\{-x)), 



(14) 

where A/' is a normalization factor. Since 6"^ ^° 

is an even function of x (since m{x) is an odd function), 
one can easily verify that 7^ = 7, i.e., the quasiparticle 
is a Majorana fermion. 

In angular momentum channels other than to, = 1, 
channel to is coupled to channel (2 — m). For these chan- 
nels, the Hamltonian is of the form 



Hr, 



(2' 



1 



mqj 



dqdq B{q,q')cl^qcl_^^, -hh.c. 



In this case there is no requirement that B{q, q') has to 
change sign upon q ^ q' since the fermion operators 
in the pairing term are from two different angular mo- 
mentum channels. Thus, generically, _B(0,0) is non-zero, 
and the spectrum is gapped. In this case the real space 
Hamiltonian resembles Eq. (2) with tp^x) the Fourier 
transform of {cmqTC2-mq)' ^^'^ mass term does not 
change sign. Hence, there is exactly one zero mode in 
the spectrum. 

For the system with a winding-number-two vortex, the 
factor e'^^ in Eq. (8) is replaced by e^'^^. This sim- 
ple modification, as we show below, gets rid of the un- 
coupled angular momentum channel. Consequently, like 
the TO ^ 1 angular momentum channels in the winding- 
number-one vortex discussed above, the real space Hamil- 
tonian is of the form in Eq. (2) where the mass term does 
not change sign. Therefore the Jackiw-Rebbi theorem 
does not apply and there is no zero energy mode in the 
spectrum. 



4 



With the modification e*^^ g^^iOn, ^^^^ function /(k + 
p) in Eq.(9) now satisfies, /(i?e(k + p)) = e2*^/(k+p). It 
follows that, 7(k + p) = e^'^^+p/dk + p|). Doing similar 
calculation as earlier we obtain, 

^ „-i{mi+m)0kf,-i{m2-m)0p J „t -i- h r' 

where 

and Vm{k,p) is the Fourier component of j^^-^- Finally, 
performing the and 6p integrals, and transforming 
m — > — m along with interchanging k and p in the third 
and fourth terms of i?2P, we arrive at the equation anal- 
ogous to Eq.(ll), 



H2P = 2Ao^ / dkdp {k^ + k'^p)^/kp Vm{k,p) 



It is clear from this equation that, in the case of the 
vortex with winding number two, no angular momentum 
channel in the pairing term decouples from the rest. 

We can generalize the above calculations to the cases 
of arbitrary odd and even winding number vortices. For 
a vortex with an odd winding number, 2n — 1, because 
of the factor e*'^^"""'^-*^^ in the pairing term, angular 
momentum channel m is coupled to channel (2n — m). 
In this case, there is always a channel, m = n, which is 
decoupled from the rest. In this channel, quite generally, 
the Hamiltonian maps on Eq. (1) with a mass term 
which is an odd function simply because of the fermion 
anticommutation relation. The index theorem we proved 
in this paper then implies the existence of a zero energy 
Majorana fermion quasiparticle. In contrast, for a vortex 
with an even winding number, 2n, angular momentum 
channel m is coupled to channel (2n + 1 — m) by the 
pairing term. Consequently, like in the case of the vortex 
with winding number two, there is no decoupled angular 
momentum channel and there is no zero mode. Note 
that in the physics discussed above the breaking of time 
reversal symmetry is crucial, since, for 2D non-chiral 
p-wave superconductors, the bulk itself is gapless. 

For s-wavc superconductors, pairing occurs between 
the fermions with opposite spins. In this case, for a vor- 
tex with winding number one, one can show that the 
pairing term couples angular momenta m and 1 — m, like 
in the uniform case discussed above. So there is no iso- 
lated channel. For a vortex with winding number two. 



the channel m = 1 decouples, and the Hamiltonian in 
this channel takes the form 

-^2 = TlT^y^ / d<l vpq c\g^cig„ + Ao / dqdq'C{q,q') 

K^T^) ^ J-A J-A 

However, in Eq. 15, because of the opposite spins of the 
fermions in the pairing term, the anticommutation rela- 
tions fail to produce an odd mass term. Thus, there is 
no zero mode at the vortices of an s-wave superconductor 
[17]. 

To conclude, we have proven an index theorem for the 
existence of zero modes in the vortex cores of chiral p- 
wave superconductors. The theorem, which is an analog 
of that in ID Dirac fermion field theory with a mass 
soliton, correctly predicts a zero energy Majorana quasi- 
particle at the core of a vortex with odd winding number, 
and no such zero mode for a vortex with even winding 
number. The formalism also correctly captures the low 
energy quasiparticle physics at the vortex cores of s-wave 
superconductors . 
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